Anisotropic quasiparticle lifetimes in Fe-based superconductors 
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We study the dynamical quasiparticle scattering by spin and charge fluctuations in Fe-based pnic- 
tides within a 5-orbital model with onsite interactions. The leading contribution to the scattering 
rate is calculated from the second-order diagrams with the polarization operator calculated in the 
random phase approximation. We find one-particle scattering rates which are highly anisotropic on 
each Fermi surface sheet due to the momentum dependence of the spin susceptibility and the multi- 
orbital composition of each Fermi pocket. This fact combined with the anisotropy of the effective 
mass, produce disparity between electrons and holes in conductivity. Hall coefficient, and Raman 
initial slope in qualitative agreement with experimental data. 

PACS numbers: 



I. INTRODUCTION 



The presence of several electronic orbitals in bands 
near the Fermi level of a metallic system provides both a 
rich set of properties and complications in revealing the 
underlying physics. Some of the most widely discussed 
examples of such systems are the recently discovered Fe- 
based superconductors with Tc up to 55K^'^ where multi- 
orbital effects cannot be disregarded. In these quasi-two- 
dimensional compounds, Fe d-orbitals form a Fermi sur- 
face (FS) consisting of nearly compensated small elec- 
tron and hole pockets. Since the sizes of the hole and 
electron FS pockets are roughly identical in the undoped 
system, one might expect a vanishingly small Hall co- 
efficient and a roughly electron-hole symmetric doping 
dependence. However, in the intensively studied 122 sys- 
tems (Ba(Fei_a;Coa;)2As2, Ba(Fei_a;Nia^)2As2) and 1111 
systems (LaFeAsOi-a^Fa; and SmFeAsOi-a^Fa^), Hall ef- 
fect measurements find that transport is dominated by 
the electrons even for the parent compounds^"—. In the 
compensated case, this result can be explained only if the 
mobilities of holes and electrons are remarkably different 
which suggests an order of magnitude disparity in relax- 
ation times, Te ^ Th'^ A similar large asymmetry of elec- 
tronic and hole scattering rates has also been suggested in 
the analysis of the electronic Raman measurements which 
can selectively probe different parts of the Brillouin zone 
(BZ) using various polarization^. Optical conductivity 
as measured by THz spectrometry provides and estimate 
of Te ~ ^Thr^ and refiectivity measurements also sug- 
gest the presence of two distinct scattering rates with a 
large disparity between themJ^— Theoretical analysis 
of the normal state resistivity p in the two-band model 
for Bai_a;Ka;Fe2As2 shows that the experimental tem- 
perature dependence p(T) can be reproduced only if one 
assumes order of magnitude larger scattering in the hole 
bandit. Finally, quantum oscillation experiments on P- 
doped systems indicate that the electron pockets have 



a longer mean free patbi^ — . It is clearly important to 
understand whether this conjectured dichotomy between 
electron and hole transport properties is real, and if it is 
universal to the Fe-based superconductors. 

There are two main sources for quasiparticle decay: 
i) electron-electron inelastic processes and ii) impurity 
scattering. We will concentrate on the first case and 
mention impurity scattering only briefiy. Experimen- 
tally, the apparent disparity in mobilities for holes and 
electrons becomes smaller as one dopes away from the 
magnetically ordered parent compound^. This suggests 
that the spin fiuctuations which also decrease upon dop- 
ing play an important role in the scattering rate asym- 
metry. Spin fiuctuations due to the nearby spin-density 
wave (SDW) state have also been considered as the most 
probable source of superconducting pairing. ^^^^ 

In this paper we study the inelastic quasiparticle scat- 
tering in Fe-based superconductors by calculating the 
scattering rate on different FS sheets within the general- 
ized spin-fluctuation theory. The self-energy is approxi- 
mated via the second-order diagrams with the polariza- 
tion operator treated in the random phase approximation 
(RPA). We show that there are two ingredients which 
provide strong anisotropy of the scattering rate. 

The most important one is that one-particle scattering 
is strongly affected by the orbital character of the ini- 
tial and final states, in analogy to orbital pair scattering 
effects which have been discussed recently^i^i leading 
to a momentum dependence of the effective interaction. 
Secondly, the polarization bubble itself is momentum de- 
pendent. The combination results in a highly anisotropic 
scattering rate on the electron Fermi surface sheets, in- 
cluding some very long lived quasiparticle states. Al- 
though our results indicate that on the average Tg is of 
the same order as r/i, the transport properties still may be 
dominated by small parts of the electron pockets where 
the lifetimes are long and the Fermi velocities are high. 
This combination causes a disparity between holes and 
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electrons in the transport properties (conductivity and 
Hall coefficient). Furthermore, analysis of the Raman re- 
sponse shows that the quasiparticle lifetime effects can be 
clearly observed in both the Big and B2g polarizations. 

A calculation of the lifetime on Fermi surface was pre- 
viously reported by Onari et al^, where the scatter- 
ing due to spin fluctuations was considered within the 
fluctuation-exchange approximation (FLEX). Our results 
reveal a similar momentum dependence of the lifetimes, 
but exhibit a much larger anisotropy due to the absence 
of self-consistency. 



II. MODEL 

We will use the 5-orbital tight-binding model of Graser 
et alS^ which is based on the ah initio band struc- 
ture calculations^^ within the local density approxima- 
tion (LDA) for the prototypical iron pnictide, LaOFeAs. 
Our interaction Hamiltonian is 

H = jj^o + ?7 riim^nimi + U' ^ ninUim 

i,m i,m<n 
^ina^lma' ^ina' Cima 

i,m<n (j,cr' 

-^J' ^ 4nt4ni^^^i^^^t' (1) 
where riim = riim^ + riimi, riima = 4ma^irna, with 

m, and a denoting site, orbital, and spin indices, re- 
spectively. The on-site intra- and inter-orbital Hub- 
bard repulsions {U and U')^ Hund's rule coupling (J), 
and the pair hopping {J') correspond to the notations 
of Kuroki et alM^ Below we will consider cases which 
obey spin-rotation invariance (SRI) through the relations 
U' = U — 2J and J' = J and those which do not. The ki- 
netic energy Hq includes the chemical potential ji and is 
described by a tight-binding model spanned by five Fe d- 
orbitals (cZxz , <^yz , d^2_y2 , d^y , d2,22_^2 The d^^ , dy^ 
and d^y bands dominate near Fermi level, as seen in Fig.[T] 
where we show the Fermi surface (FS) which arises from 
Hq in the one-Fe Brillouin zone. For the electron- and 
undoped systems the FS consists of two small hole pock- 
ets ai and a2 around the F = (0, 0) point, and two small 
electron pockets j3i and j32 around the X = (tt, 0) and 
Y = (0, tt) points, respectively. Upon hole doping a new 
hole FS pocket, 7, emerges around (7r,7r) point, which 
has been shown to strongly affect the pairing state^iiSi. 



III. METHOD 

The leading non- vanishing contribution to the quasi- 
particle scattering rate 1/r comes from the imaginary 
part of the second-order self-energy diagram (Im S) with 
the polarization bubble (see Fig. [2j). To take scattering 
from spin fluctuations into account we renormalize the 



71 


© 


dxz — 

dyz — 

Cxy 


71 


© 




x=0.03 






x=-0.08 















71 71 



FIG. 1: (Color online) Fermi surface for electron doped (dop- 
ing X = 0.03, left) and hole doped (doping x = —0.08, right) 
systems calculated within the 5-orbital model. 




FIG. 2: Orbital (a) and spin (b) structure of the second 
order diagram for the self-energy in the multi-orbital sys- 
tem, T^nn' (k, cj). Interaction lines contain four orbital indices, 
U — Unr ■ Shaded bubble denote the RPA susceptibility, 
x{q) = X^zi^i^q)- Incoming and outgoing indices n and n 
carry the same spin a. XI5 X2, and xs are the different sus- 
ceptibility channels, see Eq. JS)), a = —a. 



bubble within the random phase approximation (RPA). 
Note that second order diagrams with crossing interac- 
tion lines are not included in Fig. [2l We have chosen to 
work in this approximation to preserve consistency with 
calculations of the spin fluctuation pairing vertex^^ . The 
bubble then represents the RPA susceptibility which in 
the multi-orbital system is xliii^^^q) with w^z^v^u be- 
ing the orbital indices, and q and cUq are the momentum 
and frequency, respectively. The same susceptibility was 
calculated in Ref. |22| and was shown to produce super- 
conductivity with an Ai^ order parameter symmetry, in 
accord with several experiments^^ and other spin fluctu- 
ation calculation^i^i^i^. Here and below the orbital 
(band) indices are denoted by Latin (Greek) letters. 

Since we focus on the lifetime effects, we consider only 
Im S, neglecting the real part of the self-energy Re S. 
The renormalization of the band structure due to the 
real part of the self-energy has been discussed in some 
detail in Refs. fsolfsil and is not considered in the present 
study. We note that our calculations are based on the 
LDA band structure which already contains important 
Hartree corrections and agrees fairly well with quantum 
oscillation experimentsii^ — . 
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There are iraportant consequences of the multi-orbital 
nature of the system which deserve comment. First, the 
single-particle noninteracting Green function is diagonal 
in band space but not in orbital space. The orbital matrix 
elements a^'"^, which describe the transformation from 
one space to another are given by Ckncr = ^ct^'^di^Xai 

A 

where (ikAa is the annihilation operator for a particle 
with band index A, momentum k and energy e^. Sec- 
ondly, the interactions in Hamiltonian ([1]) have a com- 
plicated orbital structure; to compactify the expressions 
we define the local matrix interaction in orbital space, 
^n/4wai4ra2^^^^3Cin(74, which accouuts for ah the quar- 
tic terms. 

The noninteracting part of the Hamiltonian, i^o, is a 
complex matrix^^ which in general has complex eigen- 
vectors aj^'^, although the eigenvalues are real. In 
order to use a simple form of the spectral representa- 
tion of the Green function below, we choose a gauge in 
which the Hamiltonian is real by performing a unitary 
transformation Hq = (t)~^Ho4^^ where (j) is the diagonal 
matrix (j) = diag (i, i, 1, 1, 1). The interaction part of the 
Hamiltonian ([1]) must then also be rotated by (f). Having 
completed the rotation, the eigenvectors and interactions 
are now real, and after calculating the diagram in Fig. [2] 
we arrive at the multi-band extension of the standard 
zero-temperature expressions for the self-energy: 

ImE««.(k,c.) = ^ J2 f/^"f/£^'<'Vk-q (2) 

q,A w,z,u,v,r,s 

X im x7M<^- e^q) [e - e - ^)] ■ 

For simplicity, we have introduced the notation s = 
(5, CTg), where s and as are the orbital and spin index, re- 
spectively. The initial and final spins an and an' , since we 
are considering the paramagnetic state, have been kept 
equal. 

The momentum dependence of the orbital matrix el- 
ements generates an effective momentum-dependent in- 
teraction from the bare local Coulomb interactions, 

V^J{k-q) = J2ug'all\, (3) 

r 

in terms of which ([2j) may be written 

ImS]^^Kk,^) = E E KU^)VnM^) (4) 

q,A w,z,u,v 

X Im xl%{k - q, - e'j [e - e - w)] . 

The effective interaction enhances the anisotropy of the 
scattering rate, as will be demonstrated below. 

We now discuss briefly the spin structure of the dia- 
gram in Fig. [2] which is important for the calculation of 
Im S using Eq. (j2j). The susceptibility can be divided into 
charge and spin channels, and subsequently into singlet 



and triplet parts: 

^ \ Xi,2 ^ \{X')11 ± \{xTJz triplet 
I X^ = \{X')11 singlet 

where ^^^d are the charge and spin parts of the sus- 
ceptibility, respectively, and r^-cr' are Pauli spin matrices. 

For the purpose of the self-energy calculation, the in- 
teractions can be grouped into three channels. If we de- 
note the incoming spins as ai and (73, and the outgoing 
as <J2 and 0-4, the channels are: (1) cri = cr2 = (13 = (74, 
(2) ai = a2 ^ a^ = (74, (3) ai ^ a2 = a^ ^ a^. Then the 
orbital part of interactions in each channel, /7i, [/2, and 
Us, are: 

(c/i)- = o {U2)Z = u mz = -u 

(Ult'a = U'-J = U' {Ust'a = -J 

(t^i):Ho {U2)t = j' mt, = -j' 
{Ui)ii = j-u' {U2)ii = j {Us)ii = -u' 

where orbital indices a ^h. 

To combine the interactions with the susceptibility, we 
first note that due to the spin structure of the diagram, 
the interaction channels (l)-(3) decouple. Second, we see 
by inspection that channels (1) and (2) couple to xi,2, 
and channel (3) couples to xs- Thus, the self-energy will 
contain the following matrix structure 

Uxtf (X UlXlUl^U2XlU2^tfiX2tf2^U2X2U^ 

(6) 

This expression by construction resolves the spin sum- 
mation and only sums over orbital indices remain. Com- 
bining it with the calculation of x^z{^->^q) ^^r a given 
doping we use Eq. [2] to obtain Im T^nn' straightfor- 
wardly. Then we convert it to a band representation, 
Im EAA'(k, cj) = aj^'^^Im E^n' (k, a;)aj^ . For the en- 

ergy range where there are no band crossings, there is a 
unique band A corresponding to the momentum k. The 
self-energy describes the scattering of the particle with 
k back to the same momentum k, and thus back to the 
same band. A' = A. For the small energies around the 
Fermi level considered, there are no band crossings, so 
the major contribution to the scattering rate in the full 
Green function in band space, G = (Gq ^ — comes 
from diagonal, A = A', matrix elements of Im E. We de- 
note them as E^(k, cj) = Im EAA(k, cj). The momentum 
sums in Eq. [4] were performed on a 256x256 grid with 
an artificial broadening in all susceptibilities of 5 meV. 
The undoped material has completely filled d^ orbitals, 
which corresponds to ne = 6. To present our results as a 
function of doping, we define it as x = ng — 6. 

IV. SELF-ENERGY 

Because inter-band transitions are negligible in the 
range of energies considered here, the calculated scat- 
tering rate follows the Fermi liquid relation E''(k, cj) oc 
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bp' ^iPT^-^ thus, some finite frequency or temperature is 
needed for non-vanisiiing results. Here, and below, the 
quantities we report will be calculated at = 20meV 
which is equivalent to T ^ 74K at zero frequency. We 
have verified numerically that our results scale as cj^, 
and that interband transitions indeed do not contribute 
at low energies. The results below are qualitatively in- 
dependent of the exact frequency chosen, since we are 
below the range of frequencies where inter-band scatter- 
ing plays a large role. 

For several dopings and few sets of interaction param- 
eters, the calculated scattering rate along the Fermi sur- 
face is shown in Fig. [3l Here, U and J are in eV and 
were chosen to be close to the SDW-instability in the 
spin susceptibility. 

We observe that the average scattering rate increases 
monotonically with doping. Fig. |4] shows the average 
lifetime for holes and electrons on the Fermi surface, as 
well as a measure of the anisotropy which we have de- 
fined as the normalized standard deviation of the life- 
time, Ar/(r), where r/c = — l/2S''(/c, cj), scaled by the 
average. We see a clear increase in the quasiparticle life- 
time on all Fermi surface sheets as the system is electron 
doped. On the electron-doped side, the average scat- 
tering rates are essentially controlled by the degree of 
nesting. As more electrons are doped into the system, 
the hole pockets shrink and the nesting between the a 
and p> sheets deteriorates. The hole-doped systems have 
a smaller lifetime due to the presence of the 7 pocket; in 
addition to (tt, 0) scattering between a and /3 sheets, new 
phase space for scattering opens up and the average rate 
increases. Since electrons on the dxy portions of the FS 
are long-lived as will be discussed below, one expects the 
resistivity due to spin-fluctuations to increase with hole 
doping. 

Aside from the overall change in scale. Fig. |4] shows 
that the ratio of electron to hole scattering rate changes 
as one goes from hole to electron doping; electrons have a 
higher average scattering rate on the hole-doped side, and 
vice versa. Although there is already an anisotropy be- 
tween the hole and electron pockets in terms of lifetimes, 
it is not enough to cause the experimentally observed 
anisotropy, as will be discussed below. With electron 
doping, electron sheets /3i and p>2 increase in size as well 
as dxy portions. Thus the number of states with long life- 
time for electrons increases monotonically. On the other 
hand, hole pockets decrease in size and the phase space 
for scattering decreases, while for small pockets the in- 
traband scattering starts to dominate. The competition 
of these two effects lead to saturation and then to de- 
crease of lifetime for holes, indicating that the intraband 
scattering dominates. 

Next, we observe a clear anisotropy in the scattering 
rate going around the Fermi surfaces as shown in Fig. [3] 
and the inset of Fig. |4l Focusing first on the undoped 
and electron-doped systems, the ^\ sheet exhibits strong 
anisotropy between the V — X and X — M directions. 
From Fig. (TJ we observe that this is where the Fermi sur- 



face orbital composition changes from d^y to dy^ charac- 
ter. There is a strong minimum in the scattering rate 
in the d^y portions of the /3 sheets; this is due to the 
above-mentioned anisotropy of the effective interaction, 
Eq. ([3]). The orbital matrix elements tend to restrict 
scattering to be maximal for intra-orbital processes. For 
the dxy electrons, there is very little phase space to scat- 
ter compared to other orbitals, see Fig. [TJ because the 
spin fluctuation scattering intensity x(q) is peaked at 
q = (7r,0). Thus, they behave more like free electrons. 
When the system is sufficiently hole doped to create the 
(<^xy ) 7 hole pocket, (7r,0) spin fluctuations couple them 
strongly to other d^y states, causing the scattering rate 
there to increase. Throughout the doping range, dxz and 
dyz states on the a pockets scatter strongly with their 
counterparts on the /3 pockets, and vice versa. 

Finally, we discuss the interaction dependence in 
Fig. [3l The top row of panels shows a case where J = 0, 
and the middle has finite J = 0.25. As the Hund's rule 
coupling J is turned on, we observe two effects. First, 
the overall scattering rate decreases (note that the color 
scale on each plot is different). This is due to the spin- 
rotation invariance (SRI) relation U' = U — 2 J, so that 
is decreased in the middle row of panels. Although 
new scattering channels open up through J itself, this 
is more than compensated by the decrease in the inter- 
orbital scattering U\ This is confirmed by the third row 
in the figure, where J is finite but the system is non-SRI 
because U' = C/, as in the first row. Here, the scattering 
rate increases for all dopings, indicating that it is indeed 
the decrease in that is the cause of the E^' decrease in 
the 2nd row. 

Secondly, we consider the effect of J on the /3 sheet 
anisotropy for the hole-doped system. When J = 0, the 
minimum scattering rate occurs near the d^y sections of 
the Fermi surfaces for all dopings. Once J is turned on, 
the anisotropy reverses, and instead a maximum scatter- 
ing rate is found on the same sections. This reversal of 
anisotropy can be explained by the same argument as 
above. When J = 0, the intra-orbital and inter-orbital 
scattering {U and U') are the same. Thus, there is a 
strong scattering from both the d^z / dyz portions as well 
as the d^y portions of the /3 sheets to the 7 pocket (of 
dxy character). Since the dxz /dyz portions additionally 
scatter to the a sheets, a stronger scattering rate occurs 
there. When J is finite, the effective inter-orbital scat- 
tering rate U' decreases through the SRI relation. Thus, 
the scattering on the dxz / <^yz portions is decreased while 
that on the d^y sections remains the same. With suffi- 
ciently large J, the anisotropy on the P sheets is reversed. 
Note, however, that this argument depends on the exis- 
tence of the 7 pocket. When the pocket is not present, 
such as in the undoped and electron doped cases, no such 
reversal occurs, and thus the d^y states have the longest 
lifetimes for the configurations investigated. 
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FIG. 3: (Color online) Imaginary part of the self-energy E at a; = 20meV along the Fermi surface for various dopings (x = —0.14, 
0.04, and 0.13 from left to right) and for three sets of interaction parameters (in eV). All reported values are in meV. Note 
that the color scale is different for each plot. First and second row interaction parameters are SRI; third row is non-SRI. 



V. COMPARISON WITH EXPERIMENT 

A. Conductivity 

We next consider the effect of the calculated scattering 
rates on the electric conductivity. The total conductivity 
is the sum of the band conductivities, a{uj) = ^axx{(jj), 

A 



f 



(7) 



where Tk = — l/2I]^(k, cj), ki?A is the Fermi momentum 
for a particular band index A, we integrate over ky which 
is the component of momentum along the FS, ^k is the 
velocity, and A^k^A — ^/\^\<iFx\ i^ momentum- and 
band-dependent density of states (DOS) at the Fermi 
level. Note that we have approximated the transport life- 
time with the one-electron lifetime Tk, neglecting forward 



scattering corrections, as well as the distinction between 
normal and Umklapp processes. Such an approximation 
can only give the crude qualitative effect of the scattering 
from spin fluctuations on the conductivity. 

To analyze the doping-dependence of the conductivity, 
we now keep the interactions constant at values which do 
not produce an RPA instability over the range of dop- 
ings considered. We evaluate the DC conductivities at 
finite temperature by replacing l/r\^{uj) in Eq. ([7]) by 
l/^k (ttT). It is important to ask which aspects of the 
doping dependence of transport arise from purely kine- 
matic effects such as carrier density and Fermi velocity, 
which evolve with doping, and which arise from interac- 
tions. To illustrate this, we first plot in the top panel of 
Fig. [5] the separate contributions to the total conductiv- 
ity from the electron and hole sheets, with an assumed 
constant relaxation time. Here the conductivities evolve 
more or less as expected with electron doping as the vol- 
umes of hole sheets shrink and electron sheets grow. On 
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FIG. 4: (Color online) Average scattering rate for holes 
(ai, 0^2,7) and electrons {/3i,/32) at cj = 20 meV and [/ = 1.0 
eV and J = 0.25 eV. Inset: Lifetime anisotropy At/{t), 
where At ((r)) is the standard deviation (average) over the 
appropriate Fermi surface. 



the other hand, it is important that the "perfectly com- 
pensated" situation of equal kinetic conductivity of elec- 
trons and holes does not occur for the undoped case, but 
rather for x o:^ —0.05 hole doping. We have indicated in 
the figure the range of doping over which the 122 systems 
display long range magnetic order, which is not included 
in the current theory, and thus where the results are not 
directly applicable. 

By contrast, the bottom panel of Fig. [5] shows the sepa- 
rate conductivities on the hole and electron FS as a func- 
tion of doping. We immediately notice that conductivity 
for electrons grows quite strongly upon electron doping. 
Quite unlike the purely kinetic case in the top panel, the 
hole conductivity varies only weakly compared to that 
of the electrons. It is this asymmetry, due to a combi- 
nation of the kinetic effects illustrated in the top panel 
of Fig. [5] and lifetime effects calculated here, which lead 
to the rapid domination of the conductivity by electrons; 
this has led transport experiments for Co-doped Ba-122 
being interpreted in terms of a 1-band model with elec- 
trons only^'^ with some validity. The feature that greatly 
affects the doping dependence is the fact that the maxi- 
mum of the Fermi velocity is precisely where the lifetime 
is largest on the electron FS sheets, namely the d^y sec- 
tions of the f3 sheets. We also calculated conductivity and 
Hall coefficient for a case where SRI is violated {U — 1.0, 
J = 0.25, U' = [/, not shown in the figure). The re- 
sult are qualitatively similar to the case where \J = 1.0, 
J = 0. 

The conductivities obtained show a large disparity be- 
tween the hole- and electron-doped sides. It is important 
to note that what we calculate here is the spin- fluctuation 
contribution to the scattering rate, i.e. inelastic scatter- 
ing. Resistivity experiments on K-doped and Co-doped 
Bal22 show that the elastic scattering is much larger 
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FIG. 5: (Color online) Top: Conductivity for holes and elec- 
trons as a function of doping x — Ue — 6 for constant relax- 
ation rate l/r = 1 eV. Bottom: conductivity for holes and 
electrons as a function of doping x for the two sets of pa- 
rameters (in eV): U = 1.0, J = and [/ = 1.0, J = 0.25, at 
effective temperature T = 74:K. The shaded region marks the 
rough experimental SDW region in 122 systems. Solid lines 
are guides to the eye. 



in the Co-doped (e-doped) samples, presumably due to 
the fact that the Co dopants sit in the FeAs plane. This 
elastic scattering will correspondingly reduce the e-doped 
side of FiglSl and thus bring the overall scattering rate 
more in line with the experimentally observed trends. We 
have not attempted to fit experiments directly due to the 
current uncertainty in the details of the dopant scattering 
potential. 

The calculated conductivity shown in the lower panel 
of Fig. [5] was obtained for interaction parameters cho- 
sen sufficiently small to show the effect of doping while 
avoiding the RPA instability. For these parameters, the 
absolute scale of a is much larger than in experiments 
on 1111 or 122 samples. Clearly increasing the overall 
scale of the interactions will increase the scattering rates 
and decrease the conductivity. However to obtain the 
observed values of the conductivity requires approaching 
the RPA instability extremely closely. We have not at- 
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tempted to fine tune the interaction strengths, but merely 
to illustrate the possible qualitative behavior. It seems 
more likely that a more complete theory will require a 
renormalization of the susceptibility akin to that seen in 
Quantum Monte Carlo (QMC) studies of the Hubbard 
model, which indicated that the RPA form of the dy- 
namical magnetic response was qualitatively correct, but 
that the "/7" driving the instability (through the RPA 
denominator) needed to be taken independent of the U'^ 
prefactor in the effective interaction^^. A similar effect 
should occur in multi-orbital Hubbard models, such that 
the overall scales of scattering rates, and degree of prox- 
imity to the instability, should not be taken overly seri- 
ously. 
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B. Hall coefficient. 

Any disparity between the scattering rates of electrons 
and holes manifests itself in the Hall coefficient 



Rh = -cr//(cj)/cr^(cj). 



(8) 



where ani^) is the Hall conduct ivit}^!^. For a multi- 
band system, (Jh{^) = ^o-hx{^) and the expression for 

A 

the band Hall conductivity has the form 
affAH = |^ / '^kiVkVk-[Tr(Mk^)-Mk^]-VkT2(w), 

(9) 

where (Mj^ ^ ^^ = % dvk^/dkf^ is the inverse mass ten- 
sor. 

Fig. [6] shows calculated Rh as a function of doping for 
uj = 20meV (the corresponding effective temperature is 
74K). One can qualitatively understand the doping de- 
pendence of Rh by analyzing the approximate equation 
for the band Hall conductivity. 



(10) 



where 1/Rx = ±enA is the Hall coefficient for an electron 
(hole) band A, and nx is the occupation of that band. For 
the simple case of two bands (hole and electron) we have 



r)2band 
Uh 



1 ^h/^h - cflln^ 



{cTh 



(11) 



Since conductivity for the hole band (Jh oc nhTh/rrih 
and for the electron band <Je oc ngTe/me with Th^e and 
rrih^e being the corresponding lifetimes and band masses, 
j^jband jg ^ decreasing function of electron doping if 
Te ~ Th and rrie ~ rrih- This is what we see in Fig. [6] 
for the U = 1.0, J = case. On the other hand, exper- 
imental data for 1111 and 122 compounds indicate that 
i?^^* is an increasing function of electron doping (i.e., 
the magnitude \R^^^^\ decreases with increasing x) away 
from the SDW state. According to the simple analysis 
of Eq. (pT]) . this may be due to (i) Tg ^ and/or 



FIG. 6: (Color online) Doping dependence of the Hall coef- 
ficient. The theoretical calculations are for two sets of pa- 
rameters (in eV): U = 1.0, J = and [/ = 1.0, J = 0.25. 
For the first set we also show result of the multi-band ap- 
proximation for Rh from Eq. (|10p . Experimental data points 
are from (i) Ref.Hfor Ba(Fei_a:Coa:)2As2 at lOOK, (ii) Ref.l^ 
and (iii) Ref . i for SmFeAsOi_a.Fa. at 125K, and (iv) Ref . 
for BaFe2(Asi_ccPcc)2 at 150K. The shaded region tentatively 
marks the experimental SDW region. Solid lines are guides 
for the eye. 



(ii) nih ^ TUe. Note that use of Eq. ([9j) gives a dif- 
ferent result from Eq. (p!Q|) due to the mass anisotropy 
across the FS which contributes to factor (ii). Factor 
(i) starts to play a role when we consider non-zero J. 
For the case of U = 1.0 and J = 0.25, Rh{x) becomes 
slightly increasing function of x for x > (Fig. [6]). How- 
ever, it is not in quantitative agreement with experimen- 
tal data. To see whether the present approach can pro- 
vide the correct slope of Rh{x)^ we artificially increased 
scattering rate on all orbit als except d^y twice, so that 
the anisotropy between hole and electron sheets becomes 
more pronounced. The resulting doping dependence of 
the Hall coefficient is shown in Fig. [71 Now the slope of 
Rh{x) is in good agreement with experimental data. 

The fact that we underestimate the disparity between 
holes and electrons by a factor of two is not very dis- 
couraging. There are several factors not included in the 
present theory. In the interest of studying the doping 
dependence, we have kept the interactions fairly low to 
avoid the instability which occurs for relatively small 
interaction strengths on the hole-doped side. Further- 
more, we have neglected impurity scattering. In multi- 
band impurity models^^^^^, the ratio of intra- to inter- 
band scattering is taken as a parameter, and the scat- 
tering rate asymmetry between electrons and holes is 
weak. One might expect that an "orbital impurity" 
model, where an impurity introduces a local Coulomb 
potential for electrons in all (i-orbitals, might produce 
a scattering rate anisotropy in k-space due to the ma- 
trix elements a^'^^, just as in the inelastic scattering case. 
By investigating simple models similar to those consid- 
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changing function of T in Fig. [51 



FIG. 7: (Color online) Doping dependence of the Hall coeffi- 
cient for three distinct cases: (1) original calculated Rh from 
Fig. m (2) the one with the artificially increased scattering 
rate for all orbitals except for d^y , Txz.yz Txz.yz/'^^ and (3) 
Rh with added constant impurity scattering l/rimp — ImeV. 
For all cases parameters are U — l.OeV, J — 0.25eV. Experi- 
mental data points (i) are from Ref. for Ba(Fei-ccCocc)2As2 
at lOOK. The shaded region tentatively marks the experimen- 
tal SDW region. Solid lines are guides for the eye. 




40 60 80 100 120 140 160 180 
T(K) 

FIG. 8: (Color online) Calculated dependence of the Hall 
coefficient on the effective temperature for several dopings x 
with U = l.OeV, J = 0.25eV. 



C. Raman response 



ered in Ref. Is^, we have concluded that both average 
elastic scattering rate asymmetry, and elastic scatter- 
ing rate anisotropy on any given Fermi surface sheet 
are small. To address the effect of isotropic impurities 
on the Hall coefficient, we introduced a constant impu- 
rity scattering with a strength comparable to the calcu- 
lated spin- fluctuation scattering rate l/rk- Since con- 
curring scattering processes add to the self-energy, the 
scattering rate is l/r^°*^^ = l/rimp + I/tic- Substituting 
^totai gqg^ ^ Rh{x) shown in Fig. [71 

for l/rimp = ImeV. Clearly, increasing disorder leads to 
a monotonically decreasing Hall coefficient with doping 
similar to Eq. [11] with Tg ^ th- Thus dirtier samples 
will show a decrease of Rh{x) with increasing electron 
doping. 

The temperature dependence of Rh deserves addi- 
tional discussion. Recent phenomenological calculations 
of the self-energy in a two-band model for the pnic- 
tides suggest that to reproduce experimentally observed 
Rh{T) one needs to assume the non- Fermi liquid behav- 
ior of the spin susceptibility^. In particular, for large 
electron dopings, Rh{T) is almost constant but for small 
X it become an increasing function of temperature^i^. 
Here we argue that the observed temperature dependence 
can be qualitatively reproduced within our Fermi liquid 
approach. The resulting Rh{T) from our calculations is 
shown in Fig. [51 Note that the band which forms the 7 
FS pocket for X < is slightly below the Fermi level for 
small positive x. Thus at finite energy or temperature 
the scattering to that band contributes to the self-energy 
and consequently to the transport properties. That is 
the main reason why Rh{T) for x = 0.03 is a rapidly 



A momentum-sensitive probe of the scattering rate is 
provided by Raman spectroscopy. In particular, one can 
extract a scattering rate F from Raman measurements 
by considering the slope of the Raman response in the 
limit as the energy loss Q 0.^^ This quantity can be 
calculated as 



= lim7V-^ / dk-^^^^ (12) 

where 7k denotes the Raman vertex related to the inci- 
dent and scattering polarizations (see e.g. Ref. IZlh . and 
Nf is the density of states at the Fermi level. Here, 
we have taken the simplest form for the Raman ver- 
tices allowed by symmetry, namely cos{kx) — cos{ky) and 
s'm{kx) sm{ky) for the Big and B2g channels, respectively 
(note that we are using the 1 Fe unit cell conventions). 
We do not calculate the Aig response due to the diffi- 
culties involved in calculating the backflow effects 4^ In 
general, the backflow correction to the Aig channel in- 
volves the full susceptibility, not just the imaginary part. 
Although this can in principle be obtained, it is compu- 
tationally expensive. 

Fig. [9] shows the lifetimes obtained from Raman scat- 
tering according to the expression above. As discussed 
in Muschler et al.^^^^^, the Big measurements probe 
the regions of the Brillouin zone containing the electron 
sheets. The B2g measurements probe the region around 
(7r/2,7r/2), where there nominally are no Fermi surfaces. 
This causes a decrease in the overall magnitude of the B2g 
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FIG. 9: (Color online) Inverse of the Raman scattering rate 
as function of doping for [/ = 1.0 eV, J — 0.25 eV. The B2g 
curve has been scaled by 20 for visibility. 

Raman signal compared to Big , as reflected in Fig. [9l 
However, the tails of the B2g Raman vertex extend out 
to the zone edges, and thus some information can nev- 
ertheless be gleaned. On the hole doped side, the hole 
pockets are large, and the B2g vertex probes the edges of 
the hole pockets. Similarly, when the system is electron 
doped, the electron pockets grow and the B2g vertex is 
thus larger there. The numerator of Eq. [12] would give 
a symmetric doping dependence; therefore, the strong 
asymmetry is due to the lifetime effects. 

We observe that the presence of the 7 pocket has a 
large effect in the Raman response, for the same reasons 
as in the conductivity above. In particular the Big signal 
shows a large increase around zero doping. In the B2g 
channel the effect is not as strong, since there sections of 
both hole and electron sheets are probed. 

VI. CONCLUSIONS 

We have shown that the quasiparticle scattering due 
to spin-fluctuations in a multi-orbital model with local 
interactions can be signiflcantly anisotropic. Two fac- 
tors which produce this effect are the orbital matrix el- 



ements, which make interactions effectively momentum- 
dependent, and the momentum dependence of the dy- 
namic susceptibility. In the particular case of our model 
for LaOFeAs, the d^y portions of the electron FS expe- 
rience little scattering due to the small scattering phase 
space in undoped and electron-doped cases, since there 
are no d^y states on the hole sheets available for scat- 
tering. This anisotropy on the electron sheets appears 
to have profound consequences for transport in at least 
some Fe-pnictide systems. We have noted that there are 
several factors which together determine the experimen- 
tally observed disparity between holes and electrons. The 
first is the longer lifetime of the d^y states on the electron 
FS sheets. Another is the fact that the maximum of the 
Fermi velocity is precisely where the lifetime for electrons 
is largest. 

Our calculations suggest that we underestimate 
slightly the asymmetry between d^y and dxz /<^yz states 
seen in the analysis of the Hall coefficient doping- 
dependence. We have discussed and critically analyzed 
factors which can provide additional anisotropy. Finally, 
we discussed aspects of the the electronic Raman scat- 
tering rate, and showed that the lifetime effects should 
be visible in both the Big and B2g channels in varying 
amounts. 
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